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I.   INTRODUCTION 

The  subcrltical  assembly  initially  was  utilized  as  an 
experimental  unit  to  determine  if  a  chain-reacting  assembly 
were  possible.  As  reactor  technology  progressed  the  assembly 
came  to  be  used  more  and  more  as  a  test  unit  for  proposed 
reactor  designs  and  materials,  and  the  majority  of  these  test 
investigations  have  been  conducted  in  the  steady  or  equili- 
brium state.  Operation  of  the  assembly  by  disturbance  of  the 
steady  state  has  led  to  new  areas  of  investigation. 

Frequent  utilization  of  the  subcrltical  assembly  as  a  test 
unit  has  come  about  naturally.  Cost  of  a  subcrltical  assembly 
is  considerably  lower  than  that  for  a  critical  assembly  since 
special  safety  equipment  is  not  needed  (11,  p.  103).  Less 
fuel  and  moderator  are  used,  and  experiments  can  be  conducted 
easier  in  the  subcrltical  than  in  the  critical  assembly. 

Disturbance  of  the  steady  state  may  be  accomplished  by 
two  methods,  reactivity  forcing  and  source  forcing.  Re- 
activity forcing  can  be  represented  as  a  change  in  the  neutron 
density  brought  about  by  a  small  harmonic  variation  In  the 
reactivity,  whereas  source  forcing  is  produced  by  a  variation 
in  the  rate  at  which  neutrons  are  being  added  to  the  assembly. 

Thus  far  the  investigation  of  the  dynamic  response  of  a 
system  has  been  limited  to  servomechanism  theory,  which 
treats  the  system  as  a  "black  box" .  Phase  and  amplitude 


response  of  the  system  in  relation  to  a  forcing  function  must 
be  determined  in  order  to  apply  servomechanism  theory.  The 
above  relationship  called  the  transfer  function  is  derived  by 
noting  the  responses  as  indicated  brought  about  by  either 
reactivity  or  source  forcing. 

It  is  believed  however  that  neutron  diffusion  in  a  sub* 
critical  assembly  will  not  produce  identical  phase  and 
amplitude  responses  throughout.  These  variations  should  depend 
upon  the  method  of  varying  the  rate  of  neutrons  from  the 
source  and  the  position  of  neutron-flux  measurement. 

The  purpose  of  this  investigation  was  threefold.  First, 
an  experimental  verification  of  the  theoretical  transfer 
function  for  a  subcrltical  assembly  was  attempted  using  a 
varying  rate  of  neutrons  from  the  source.  Secondly,  neutron 
diffusion  in  the  subcrltical  assembly  that  shows  variance 
from  the  predicted  identical  phase  and  amplitude  response 
throughout  the  assembly  was  investigated.  The  last  area  of 
investigation  was  an  attempt  to  determine  assembly  character- 
istics by  using  a  sinusoidally  varying  source. 


II.  REVIEW  OP  THE  LITERATUBE 

A  reasonable  amount  of  literature  la  available  on  the 
development  and  usage  of  the  transfer  function  of  a  critical 
assembly.  Original  work  In  the  field  was  carried  out  by 
Owens*  Crever  and  Pigott  (12)  In  19*9,   In  which  an  automatic 
control  unit  was  used  to  compensate  for  reactivity  variations 
In  order  to  maintain  the  desired  power  level  In  a  reactor. 
Their  work  showed  that  the  assembly  transfer  function  ex- 
hibited proper  characteristics  to  describe  the  reactivity 
variations  in  the  assembly. 

Later  that  same  year,  Franz  (4)  developed  the  general 
form  of  the  transfer  function  of  a  nuclear  reactor  by  apply- 
ing the  standard  techniques  of  servo  theory.  He  applied  the 
transfer  function  to  an  electronic  assembly  simulator  to  pre- 
dict assembly  response  to  a  step-function  change  of 
reactivity. 

The  practical  verification  of  using  a  subcrltlcal  assembly 
In  testing  reactor  components  was  shown  by  Axtmann,  Dessauer, 
and  Parkinson  (1)  In  1)55.     A  chain-reacting  'test  assembly" 
was  utilized  to  test  materials  by  noting  effects  on  the  kinet- 
ic behavior  of  the  assembly.  Identical  tests  were  conducted 
in  a  subcrltlcal  assembly,  and  the  results  were  nearly 
identical  with  the  subcrltlcal  assembly  producing  results 
faster,  cheaper  and  safer* 


Qlasstone  (})  briefly  discussed  the  transfer  function  in 
a  subcritical  assembly.  In  his  development  he  assumed  a 
sinusoidal  variation  in  source  strength  and  used  electrical 
analogs  of  reactor  quantities  to  indicate  the  relationship  of 
electrical  parameters,  with  the  reactor  parameters,  reactivity 
forcing  and  source  forcing. 

Weinberg  and  Schweinler  (J&)  developed  the  equations  for 
the  response  of  a  critical  chain-reacting  assembly  to  a 
thermal  neutron  absorber  which  is  oscillated  back  and  forth 
inside  the  assembly.  They  further  stated  that  at  frequencies 
which  are  low  compared  to  the  reciprocal  periods  of  the  de- 
layed neutrons,  the  neutron  density  in  the  assembly  rises 
and  falls  as  a  whole.  However,  if  the  frequency  is  much 
greater  than  the  reciprocal  period  of  the  delayed  neutrons, 
a  spherical  neutron  wave  emanates  from  the  vicinity  of  the 
oscillator  and  is  attenuated  in  an  exponential  manner  with 
distance  as  it  is  propagated. 

Weinberg  and  Wigner  (16)  extended  this  theory  to  sub- 
critical  assemblies.  It  was  stated  that  the  neutron-wave 
wave  length  is  related  to  the  material  buckling.  They  sug- 
gested that  the  exponential-experiment  method  of  determining 
material  buckling  could  be  performed  by  measuring  the  proper- 
ties of  neutron  waves  established  by  a  localised  oscillating 
neutron  source.  A  particularly  desirable  feature 
of  this  method  is  that  an  experiment  to  determine  material 


buckling  can  be  performed  In  an  assembly  of  smaller  dimensions 
than  that  required  by  the  normal  exponential  experiment. 
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III.   LIST  OF  SYPBOLS 

A    Activity  counts/min. 

a    Width  of  assembly  in. 

a.l.  Attenuation  length  in.,  cm. 

Bz      Buckling  cm.~a 

@       fraction  of  average  group  of  delayed 
neutrons 

^  Fraction  of  ifcn  group  of  delayed  neutrons 

C    Concentration  of  average  group  of  delayed 
neutron  precursors 

Ce   End-KJorrection  factor 

Cft   Harmonic  -correction  factor 

ft   Concentration  of  l">  group  of  delayed 

neutron  precursors 

D    Thermal  neutron  diffusion  coefficient       cm. 

f    Thermal  utilisation 

Y        Inverse  relaxation  length  in."1 

K.    Cons  bant 

k        Effective  jsultiplication  factor 

fcc^     Multiplication  factor  for  an  infinite 

assembly 

L        Neutron  diffusion  length  in  assembly  cm. 

X        Average  neutron  lifetime  in  a  finite 

assembly  sec  • 

$        Vteoxk  neutron  lifetime  in  an  infinite 

assembly  sec . 


s 


X  Decay  constant  of  average  group  of  de- 
layed neutron  precursors 

X±  Decay  constant  of  the  i   group  of  de- 
layed neutron  precursors 

ii-  Migration  area 

n  Number  of  neutrons  of  thermal  energies 

(J  Frequency  of  oscillation 

p  Resonance  escape  probability 

$  Thermal  neutron  flux 

X  Macroscopic  cross  section 

t  Time 

T  Fermi  age 

v  Thermal  neutron  mean  velocity 

z  Vertical  distance  above  base  of  assembly 
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IV.   THEORETICAL  ANALYSIS 

A.  Transfer  Function 

The  transfer  function  is  defined  as  the  ratio  of  the 
phase  and  amplitude  response  to  a  sinusoidally -varying  forcing 
function.  This  response  for  a  subcritical  assembly  can  be 
obtained  from  the  time -dependent  general  diffusion  equation. 
Thermal  neutrons  only  will  be  considered  since  the  experi- 
mental equipment  is  designed  to  vary  the  rate  at  which  thermal 
neutrons  enter  the  assembly.  The  basic  equation  is  (6, 
p.  101) 

J  V  ^  /i  -  Za  p  •*•  S  »  — r~  m ~—     Eq.  1 

where  0   *  n  v,  and  v  is  assumed  to  be  constant. 

The  source  term  can  be  separated  into  three  parts.  The 
first  term  is  due  to  the  prompt  neutrons  resulting  from 
fission.  Fermi -age  theory  yields  this  part  of  the  source 

term  to  be  (^,  p.  226) 

Sp  -  (1  -  0)  k^  Za  0  e  '  Eq.  2 

where  (3  is  the  fraction  of  all  the  delayed  neutrons. 

The  second  part  of  the  3ource  term  is  due  to  the  delayed 
neutrons  resulting  from  decay  of  various  fission  products. 
A  total  of  six  groups  of  delayed  fission  neutrons  are 
utilized  in  expressing  this  part  of  the  source  term.  Each 
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group  has  its  own  value  of  X^,    Pj  and  Cj,.  The  net  rate  of 
formation  of  the  delayed -neutron  precursors  Is  then 


di 


d 


(3i  ~-    Za  pf  -  Ai  Ct  .      Eq.  3 


The  delayed -neutron  source  term  using  Fermi -age  theory  can  be 
expressed  as  (5*  p.  227) 

S^n  «  P  •     ^/_  ^1  ci  E<** 

1*1 

where  21    Ai  C^  is  the  decay  rate  of  the  precursor  which  Is 

i*»l 

numerically  equal  to  the  rate  of  production  of  the  delayed 
neutrons.  The  Feral-age  for  the  delayed  neutrons  has  been 
taken  equal  to  that  of  the  prompt  neutrons  J$j  p.  22?). 

The   last  part  of  the  source  tens  Is  due  to  the  external 
source  which  in  this  analysis  supplies  the  source  forcing 
needed.  It  is  written  as 

Se  •  $0  +  sA  ©^  ^  Kg*  g 

With  substitution  of  the  above  expressions  of  the  source 
term  into  Equation  1  the  following  equation  is  obtained 

»  v*  p  -  ra  *  +  u-  e)  koo  za  *  e"B*r 

♦  p  e  ***r  £  ^c^So+s^^^ 

i*i 

*  .J diL  •  Eq-  6 

v   at 
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The  space  and  time  variables  arc  considered  separable 
when  k  remains  constant  (5*  p.  227).  This  yield3  an  ordinary 
differential  equation  for  Equation  o.  If  the  subcritical 
system  Is  large,  then  va  0  «  -Ba  0   (6,  p.  361). 

Let  /0  «=  Eq.  i 

and  L2  «=  D/  £a  Eq. 

and  divide  Equation  6  through  by  £a,  to  obtain  the  result 

-  (LaB2  +  I)  0   +  (1  -  (3)  k^  e~B  T  0 


S     ej  ut                    -  * 

y                      *•     d  tP       < 
^-a 

Eq.  9 


With  division  of  Equation  9  through  by  (1  +  L2B2)  and  the  use 

of  the  relationships 

-BaT 

k  »  — I! _ Eq.  10 

1  +  L2B2 


and 

Eq.  11 


i. 


1  +  L^E2 
Equation  9  becomes 


T 
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k     6 


♦  (l  -  (3)  ic  V  +  -f £ Z  At  Ci 


k°° 


So^S,^^     -  d  ft 


Za  (1  +  L*B*)      *  t 


Eq.  12 


With  division  K    JL  \     *&   utilization  of  #  *  n  v,  Equation 
12  Is  reduced  to 

n  v     p    k 
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Ck(i  -  0)  -  13  —a—  +  — = —  -TT- —  z_  At  ^ 
SQ  ♦  S.  e^  "  *      d  n 


v  H  .   .  Eq*  13 


Zai(l  ♦  L*Ba)      d  t 
With  rearrangement  Equations  IC  and  11  become 


T1 

M  ««__  i|    2«  v  Eq*  14 


— *- 


^     CO 


and 


X 
A  -  i(l  *  L®B*)  « 


or 


1/v  «  Z  ^  (l  +  LaBa)   .  Eq.  15 


Equation  13  with  use  of  the  above  two  substitutions  and 
division  through  by  v  yields 

[k(l  *  (3)  -  1)  —  +  p  e~^T  Y     A|  C± 

+  S0  ♦  SA  ej  u  Z   -  ~— -  .  Eq. 
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In  a  similar  manner,  Equation  3  with  the  use  of  Equations 
10  and  13  is  found  to  be 


d  C«         k  n 


d  t 


3«r"~  "  ^i  Ci  •   ;   Sq.  17 


H 


Solving  for  X  ^  C^  and  substituting  into  Equation  16  yields 


[k(l  -(3)   -  13  -j-  +  P  e~B*T    £ 


n 


1-1 


(3<   k  n 


p  e 


:r~'f 


d  C± 
d  t 


•4-   Sq    +    Si    @ 


j  CJt 

d  t 


Eq.   18 


or 


-  p  e 


-B2T 


y  d  cj 

i=i   d  t 


■f  SQ  +  Sx  e 


J  6J  t 


Eq.   19 


since 


6       (3      kn  (3kn 

/  '  ..  SB  '  • 

1»1  X  A 


Solving  thi3  for  n  gives 


n  » 


(k-1) 


P  e 


-B*T 


d  ft 


d  n 


-*.  ^-X 


d  t 


d  t 


l  (SQ  +  Sx  eJ  ^  *  ) 
U-k) 


Bq.   20 
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In  thi3  particular  case  the  value  of  k  will  remain 
constant.  Neither  assembly  geometry  nor  fuel-loading  changes 
will  take  place  in  the  experiment.  However  both  the  neutron 
density  and  precursor  concentration  will  vary  as 

n  m  n0  +  nL   e^  ^  *  Eq.  21 


and 


ci  *  coi  ♦  cii  eJ  U  t  E<2-  22 


where  nx   and  Cx   are  the  complex  amplitudes.  With  the  aid  of 
Equations  17 ,   21  and  22,  Equation  20  yields 


,JWt  , 


n0  +  nx  e*    -  (k>1) 


-  w  T  V 
p  ©    2. 
i«i 


— r —  (no  +  nx  e     ) 

t  p  •*  T 

-  A|  (Col  ♦  ***  Jut)l    ^^~(n0  +n,  t*^*) 

i  (S0  tS.e^M 

♦ .  Eq.  23 

(1-k) 

The  complete  solution  to  Equation  23  involves  using  the 
six  groups  of  delayed  neutrons  each  with  its  discrete  value 
o£  X  i#  (3i»  and  Ci.  The  solution  using  an  average  group  of 
delayed  neutrons  will  be  shown. 

Equation  23  can  be  separated  into  steady  and  sinusoidal 
components,  and  each  constitutes  a  valid  equality.  The 


;: 


u 


transfer  function  uses  only  the  sinusoidal  component ,  which 
with  group  notation  eliminated,  la 


n,  a 


6o> 


i 


(3 


di-v  %        /.it 


j  6J  t 


(1-k) 


Eq.  2^ 


Again  usliu  Equation  17  and  substituting  Equations  21 
and  22  gives 

k  (3 


jrp-  (no  ♦  »i  «J  w  *)  -  *  <co  +ck**u  ») 


♦-2*-(Co  *c4  •^^t)   ■ 


The  sinusoidal  portion  or  which  becoiaes 


Eq.  25 


k  0 


•B*r 


nx  a*1  ^  r  «     X0k  eJ  w  ' 


*   4     6J    Cx    ©^    U  l  , 


£q.   26 


which  reduces  to 


(3  k  n 


U     +    jW)     C; 


Pi§ 


8  r 


Eq*  27 


Substituting  Equation  27  into  Equation  24  yield* 


lb 


nt  e 


u) 


(3  k   nx  e«*  ^  * 
(k-1) 


(k-1) 


(3  k  tti 


j?St  e«*  u  t 
(1-k) 


i  2unx   e 


jut 


(k-1) 


Eq.  26 


By  definition,  the  amplitude  portion  of  the  transfer 
function  is  the  ratio  of  the  response  to  the  sinusoidally  - 
varying  driving  source.  The  ratio  can  be  expressed  as 


n 


S 


i*k 


Eq.  29 


or  from  Equation  28 


[1  + 


&  k 


(1-k) 


A   k    (3 


j?  4  cj 


Eq.   30 


(  A  tjw)(l-k) 


(1-k) 


■] 


Rearranging  Equation  30  by  putting  the  right  side  denominator 
under  the  common  factor  ( X  +  j£o)  (1-k),  multiplying  through 
fcy  A  ♦  J6J#  then  dividing  by  A  and  simplifying  gives  the 
desired  form 
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n. 


Sx 


1-lc 


A(l-k) 


1  + 


JCJ 


4cj  r, 


(Mc) 


Eq.  31 


B.  Attenuation  Length  of  Neutron  Wave 


When  the  rate  at  which  neutrons  entering  a  subcritlcal 
assembly  is  oscillated  in  a  sinusoidal  manner  at  a  frequency 
considerably  greater  than  the  reciprocal  period  of  the  de- 
layed neutrons  a  neutron  wave  is  propagated.  With  an  average 
delayed  neutron  period  of  approximately  10  seconds,  fre- 
quencies greater  than  one  radian/second  should  produce  these 
neutron  waves*  Frequencies  considerably  less  than  the 
delayed -neutron  reciprocal  period  allow  the  delayed  neutrons 
to  remain  in  equilibrium  with  the  prompt  neutrons  (14, 
p.  100).  However,  at  higher  frequencies  the  assembly  be- 
haves as  if  there  were  no  delayed  neutrons  and  the  delayed 
neutrons  merely  contribute  to  the  background.  Hence  the 
neutron  wave  has  characteristics  that  are  dependent  upon 
the  "prompt"  response  of  the  assembly.  These  characteristics, 
wave  length,  velocity,  and  attenuation  length  are  all 
dependent  upon  the  frequency  of  oscillation. 

Weinberg  and  Wigner  (16,  p.  ^37)  developed  equations 
for  the  three  neutron-wave  characteristics.  lUe  equation 


17 

for  the  characteristic  that  is  of  Importance  here,  attenua- 
tion length,  is 

a.l.  * 

J[(kp  -  l)a  +  icjj0)*}*/*  -  (top  -  D]  */« 

where  cj  is  greater  than  one  radian/second.  Attenuation  length 
is  defined  as  the  distance  for  wave  amplitude  to  decrease  by 
a  factor  of  ©. 

If  a  neutron  wave  were  propagated  by  using  a  high- 
frequency  sinusoidal  source  oscillation  then  attenuation 
length  could  be  determined*  Amplitudes  of  the  sine  wave 
in  selected  vertical  positions  could  be  measured,  and  a  value 
for  attenuation  length  would  then  be  determined  from  the 
equation  A  ■  K  e"*'a   ,  where  A  is  the  amplitude  at  z,  and 
K  is  a  constant.  A  semi-log  plot  of  amplitude  versus  verti- 
cal position  will  yield  a  slope  of  1/a.l.  This  value  of 
1/a.l.  is  the  "inverse  relaxation  length"  of  the  complex 
buckling  for  the  frequency  investigated.  If  the  phase  angle 
is  small  at  a  given  frequency,  the  attenuation  length  at  that 
frequency  would  be  approximately  equal  to  relaxation  length 
used  to  calculate  material  buckling.  Therefore,  the  complex 
buckling  determined  by  this  value  of  attenuation  length 
would  be  approximately  equal  to  the  material  buckling. 


i.. 

V.  EXPERIMENTAL  EQUIPHEHT 
A.  Subcrltieal  Assembly 

The  subcrltieal  assembly  used  for  the  experimental  in- 
vestigation is  shown  In  Figure  1  with  the  east  cover  removed. 
Cylindrical  rods  of  graphite  measuring  7  Inches  in  diameter 
were  machined  to  a  square  from  6  inches  across,  leaving 
rounded  corners.  These  blocks  of  AGR-grade  graphite  were 
then  stacked  In  ten  columns  nine  rows  high  to  form  the  lower 
section  of  the  assembly.  For  the  top  five  rows  cylindrical 
graphite  rods  6-3A  Inches  in  diameter  were  machined  to  a  5 
inch  by  6  inch  cross  section  with  rounded  corners.  The 
space  available  between  the  corners  of  adjoining  blocks 
was  utilised  for  insertion  of  fuel  elements  or  measuring 
equipment . 

The  graphite  blocks  were  assembled  upon  a  plywood  sheet 
supported  by  a  wooden  pedestal.  As  shown  in  Figure  1,  the 
pedestal  was  divided  into  three  accessible  spaces  approxi- 
mately one  foot  high.  The  two  outermost  spaces  contained 
open-topped  aluminum  tanks ,  filled  with  water,  which  ex- 
tended the  width  of  the  assembly.  They  served  to  moderate 
the  neutrons  and  to  reduce  radiations  coming  directly  from 
the  sources  located  in  the  oscillator  unit  beneath  the 
assembly.  Figure  2   shows  the  oscillator  unit  and  selective- 
count  device  which  was  placed  in  the  center  access  section 
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Figure  X,     The  suberitlcal  asasiably 
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Figure  2.  The  oscillator  unit  and  selective- 
count  device 
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beneath  the  assembly. 

All  four  sides  and  top  of  the  assembled  graphite  were 
covered  with  10-mil  cadmium  sandwiched  between  a  0. 375-inch 
thickness  of  plywood  and  0.125  Inches  of  masonlte.  Approxi- 
mately 95  percent  of  all  thermal  neutrons  were  stopped  by 
this  cadmium  sheet,  giving  a  "black  boundary"  to  thermal 
neutrons.  A  wooden  framework  braced  the  north  and  south 
sides  of  the  assembly  and  both  the  east  and  west  faces  were 
removable  for  easy  access. 

Natural  uranium  in  the  form  of  cylindrical  rods  1.0 
Inches  in  diameter  and  o.O  inches  long  served  as  the  fuel. 
These  rods  were  encased  by  2S  aluminum  cans  with  a  0.010- 
inch  wall  thickness  and  end  caps  C.200  inches  thick.  Each 
can  of  uranium  was  then  helically  wound  with  a  spacing  wire 
and  placed  into  a  61S  aluminum  tube,  62  inches  long  with  an 
outside  diameter  of  1.375  inches  and  a  wall  thickness  of 
0.035  inches.  Approximately  10  feet  of  the  2S-aluoinum 
spacing  wire  was  needed  for  each  tube.  A  total  of  seven 
canned  uranium  rods  were  placed  in  each  tube  to  constitute 
one  fuel  element.  With  the  fuel  elements  placed  in  every 
other  hole  as  shown  in  Figure  1  an  5.5-inch  square  lattice 
in  the  lower  region  of  the  assembly  was  formed.  This  lattice 
configuration  was  maintained  throughout  the  experimental 
investigation. 
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B.  Oscillator  Unit 

The  oscillator  unit  and  select,   -.aunt  device  Is  shown 
In  Figure  2.  The  drum-type  assembly  at  the  right  end  is 
the  oscillator  unit  itself  and  la  show*  s^        ly  in 
Figure  3.  This  unit  was  designed  to  contain  the  five  one- 
curie  Pu-Be  sources  and  to  vary  the  rate  at  which  thermal 
neutrons  were  emitted  from  the  source  in  a  sinusoidal 
manner. 

Each  cylindrical  source  measuring  1  inch  in  diameter 
and  l-3/o  inches  high  was  placed  in  the  source  retainer. 
A  non-retating-paraffin  jacket  surrounding  the  retainer  was 
constructed  of  1/16-lnch  aluminum  sheet  rolled  into  a 
diameter  of  7-3/S  inches.  2nd  caps  of  1/6-inch  aluminum 
were  attached,  giving  an  overall  length  of  **•§  Inche 
Twenty-five  pounds  of  white  paraffin  were  poured  into  the 
Jacket  through  a  filler  cap,  and  precautions  were  taken  to 
reduce  the  possibility  of  shrinkage  voids,  Moderation  to 
supply  the  thermal  neutrons  desired  was  provided  by  the 
three  Inches  of  paraffin.  The  paraffin  jacket  and  the  hear- 
ings for  the  pattern  cylinder  were  both  supported  by  the 
source  retainer  which  extended  to  the  support  brackets. 

An  outer  rotating  cylinder,  the  pattern  cylinder,  was 
also  f erased  from  1/8-inch  aluminum  sheet.  Its  dimensions 
were  an  outside  diameter  of  8  inches  and  a  length  of  15 
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ROTATING     CYLINDER 
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JACKET 
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(360°) 


Figure  3.  Oscillator  unit  and  cadmium-pattern  detail 
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Inches.  End  plates  made  from  1/4-inch  aluminum  were  fitted 
and  serened  to  a  collar  fitted  on  Ita  bearing  outer-run 
casing.  The  inner -run  casing  of  the  bearing  rested  on  the 
source  retainer  for  support.  To  give  easy  access  to  the 
inner  section  of  the  pattern  cylinder  the  end  plates  were 
attached  by  screws. 

A  sine  pattern  with  the  dimensions  as  shown  in  Figure  3 
was  cut  from  10-mil  cadmium.  Three  thicknesses,  giving  a 
total  of  0,030  inches,  were  cemented  to  the  inner  surface  of 
the  pattern  cylinder.  !4ost  of  the  thermal  neutrons  were 
stopped  and  a  small  fraction  of  the  fast  attenuated  by  this 
thickness  of  cadmium  (10,  p.  12). 

Point  A  of  the  pattern  was  the  wide-open  position,  allow* 
ing  the  maximum  number  of  thermal  neutrons  to  escape  from  the 
oscillator  unit.  The  closed  position,  point  B,  stopped  most 
of  the  thermal  neutrons.  When  point  A  was  at  top-dead 
center,  the  oscillator  unit  was  at  the  zero-degree  position. 
The  Is'ft -degree  position  occurred  when  point  B  was  at  top-dead 
center.  Hence  the  exact  position  of  the  sine  pattern  of 
the  oscillator  unit  ma  designated  by  a  given  angle  from  0 
to  360  oscillator  degrees. 

Measurements  taken  within  1.3  inches  of  the  oscillator 
unit  at  the  two  extreme  positions  described  above  showed 
that  the  paraffin  jacket  moderated  a  satisfactory  percentage 
of  neutrons  to  be  used  in  the  experimental  investigation. 
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The  value  obtained  at  position  A  minus  that  at  B  was  56  per- 
cent of  the  value  at  position  A. 

c.  Counting  Apparatus 

In  order  to  measure  the  neutron  flux  during  a  designated 
number  of  oscillator-unit  degrees  a  selective -count  device 
was  constructed.  A  positive ,  non-slip  linkage  had  to  be 
supplied  from  the  oscillator  unit  to  a  position  removed  from 
the  neutron  sources.  The  equipment  needed  to  select  the 
number  of  oscillator  degrees  in  which  counting  would  take 
place  was  located  at  this  "removed"  position. 

A  bicycle  sprocket  was  bolted  to  the  end  plate  of  the 
pattern  cylinder.  A  110-inch  bicycle  chain  was  used  to 
connect  this  sprocket  on  the  oscillator  unit  to  a  sprocket 
of  the  same  diameter  mounted  on  a  countershaft  at  the  removed 
position.  Position  by  degrees  of  the  oscillator  unit  was 
then  transmitted  unchanged  to  the  "removed"  position  where 
work  could  be  performed  without  the  intense  radiation. 

The  countershaft  was  a  polishing -head  spindle  which  was 
mounted  on  a  3/^-inch  plywood  extension  connecting  the 
oscillator  unit  to  the  remainder  of  the  equipment.  At  one 
end  of  the  spindle  shaft  a  6-inch  V-belt  pulley  was  locked 
to  the  shaft,  and  the  bicycle  sprocket  was  in  turn  bolted 
to  this  pulley.  At  the  opposite  end  of  the  shaft  was  locked 
a  plywood  plate,  6-3/^f  inches  in  diameter.  The  polishing- 
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head  spindle  was  driven  by  a  1/4-inch  Black  and  Decker 
electric  drill  connected  to  it  by  a  leather  belt.  Rotation 
frequency  was  controlled  by  varying  the  voltage  to  the  drill. 
In  order  to  improve  the  frequency  control  at  alow  apeeda 
(below  £G  HPM)  a  second  polishing -head  spindle  was  interposed 
as  shown  in  Figure  2.  'These  two  configurations  then  gave  a 
speed  range  for  the  oscillator  unit  of  from  k  to  1000  BP&. 

A  1/4-lnch  plywood  plate  of  6.14-inches  diameter  with  a 
30-degree  cam  hump  was  fabricated.  To  prevent  wear  a  metal 
strip  was  cemented  to  the  cam.  A  similar  10-degree  cam  was 
also  constructed.  The  cam,  with  a  1/2-inch  centering  hole, 
was  held  securely  by  a  nut  against  the  plywood  plate  on  the 
spindle  shaft.  A  compass  card  oriented  to  the  exact 
oscillator  degree  was  secured  to  this  plywood  plate.  The 
cam  was  set  at  any  desired  oscillator  position  by  loosening 
the  cam-holding  nut. 

As  shown  In  Figure  2   a  microswltch  was  positioned  adja- 
cent to  the  cam.  The  switch  was  wired  in  the  normally-open 
position  to  the  preset-time  relay  of  a  Radiation  Instrument 
Development  Laboratories,  model  206,  scaler.  The  cam  was 
positioned  by  use  of  the  compass  card  to  make  contact  with 
the  switch  at  a  specified  oscillator  degree.  When  the  sine 
pattern  of  the  oscillator  unit  reached  this  degree  the  cam 
made  contact  and  the  switch  was  closed.  This  closed  the 
preset-time  relay  allowing  the  scaler  to  begin  to  count. 
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When  the  mlcroswltch  broke  contact  the  relay  of  the  scaler 
opened,  ceasing  the  count. 

A  Radiation  Counter  Laboratories  boron  trif loride  pro- 
portional counter 9  model  10pQ3*  was  inserted  midway  into  the 
assembly  through  holes  cut  into  the  oast  cover.  It  in  turn 
was  connected  to  the  Radiation  Instrument  Development 
Laboratories  scaler. 

The  oscillator  unit  and  selective-count  device  was  in- 
serted from  the  south  end  of  the  assembly  until  the  oscillator 
unit  was  at  the  center.  As  shown  in  Figure  k   this  placed 
only  the  oscillator  unit  Itself  beneath  the  assembly}  the 
remaining  equipment  was  outside  to  facilitate  changing  the 
earn  setting. 


28* 


Figure  *u  Experimental  set  up 
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A*  Description  of  Typical  Hun 

The  fuel  elements  were  loaded  into  the  assembly  as  shown 
In  Figure  1*  Measurements  were  taken  only  in  the  center  of 
the  assembly-  Holes  drilled  on  the  east-cover  centcrline 
matched  assembly  holes  between  the  graphite  blocks.  The 
holes  were  one  foot,  two  feet,  three  feet,  four  feet  ana 
four  feet  eleven  Inches  above  the  base.  Hereafter  these 
holes  will  be  referred  to  as  holes  one  through  five 
respectively. 

She  neutron  detector  was  inserted  midway  Into  the  assembly 
and  the  oscillator  unit  set  at  0  degrees.  With  the  30- 
degre©  cam  set  at  a  mean  position  of  0  degrees,  initial 
contact  with  the  mieroswitch  was  made  at  3^5  degrees. 

Oscillator  speed  was  adjusted  to  the  desired  HFH  by  re- 
ferring to  a  stop  watch  and  counting  revolutions.  When  de- 
sired speed  was  obtained  the  watch  was  started  and  the 
preset-time  switch  on  the  scaler  was  actuated.  Huns  varied 
In  duration  from  one  to  five  minutes  depending  on  position 
and  oscillator  speed,  A  run  was  completed  when  a  designated 
number  of  counting  periods  had  elapsed  rather  than  ending  at 
an  exact  time  interval.  For  example,  a  run  for  3  minutes 
at  18  R.m  requires  $b  counting  periods j  when  $4  periods  had 
elapsed  the  preset-count  switch  was  turned  off  even  though 
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the  time  may  have  been  a  few  seconds  on  either  aide  of  3 
minutes.  It  was  difficult  to  control  the  speed  exactly  and 
this  procedure  eliminated  varying  numbers  of  counting  periods 
in  a  single  series  of  runs* 

The  counts  were  recorded  from  the  scaler j  the  drill 
turned  off;  the  cam  reset  to  a  mean  position  of  30  degrees 
and  the  entire  counting  procedure  repeated.  The  counting 
procedure  continued  in  30-degree  increments  until  the  entire 
360  degrees  had  been  surveyed  yielding  one  single  series  of 
runs. 

Speeds  (frequency  of  oscillation)  used  in  the  investiga- 
tion were  selected  at  regular  intervals  on  a  logarithmic 
scale.  Runs  were  attempted  throughout  the  entire  range  of 
the  equipment.  The  lowest  constant  speed  attainable  was 
4  RPM.  An  upper  limit  on  speed  was  imposed  by  the  cam  and 
microswitch  arrangement,  for  in  speeds  beyond  160  EfH  the 
switch  hammered  the  cam  as  the  switch  spring  forced  the  idler 
back;  to  the  surface  of  the  cam.  As  runs  were  made  in  the 
higher  hole  positions  at  high  speeds  the  data  became  unusable. 
This  was  caused  by  the  exponential  attenuation  of  neutron 
flux  up  the  assembly  and  a  consequently  smaller  difference 
from  the  highest  to  the  lowest  count  rate  in  a  series  of 
runs.  When  this  occurred  the  counting  rates  when  plotted 
•gainst  degrees  did  not  yield  a  satisfactory  distribution. 

Original  data  for  the  various  runs  completed  is  shown  in 
Tables  5  through  13  in  the  Appendix. 
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B.  ife termination  of  Phase  Angle 

The  count  rates  for  a  single  series  of  runs  were  plotted 
on  graph  paper.  Standard  deviation  for  each  point  was  also 
plotted  to  facilitate  drawing  a  representative  curve  through 
the  points.  All  the  plots  that  showed  a  trend  in  distribu- 
tion gave  a  sinusoidal  result*  Figure  5  shows  a  plot  of  the 
data  for  all  holes  at  2. 51  radians/second,  and  as  mentioned 
previously  some  unusable  data  was  taken  and  the  basis  for 
declaring  it  unusable  was  large  scatter  in  the  plotted 
curves  * 

From  the  general  trend  of  the  plot  and  especially  from 
the  position  of  the  highest  and  lowest  counting  rates  the 
phase  angle  for  that  series  of  runs  was  determined.  To  ensure 
that  an  accurate  phase  angle  was  be Ins  obtained  a  10-degree 
earn  was  utilised*  Five-minute  counts  were  taken  every  ten 
degrees  in  the  region  of  the  maximum  counting  rate,  however 
use  of  the  10-degree  cam  was  restricted  to  the  lower  speeds. 

It  was  assumed  in  determining  the  phase  angle  that  the 
oscillator  unit  when  in  the  O-degree  position  during  rota- 
tion and  the  counting  equipment  would  yield  the  highest 
counting  rate,  to  validate  this  assumption  an  entire  series 
of  calibration  runs  with  the  neutron  detector  held  1*5 
Inches  fro*  the  oscillator  unit  In  position  beneath  the 
assembly  was  made.  Table  13  In  the  Appendix  and  Figure  5 
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indicate  the  results  which  show  that  there  was  no  "built  in" 
phase  angle. 

C.  Determination  of  2f  by  Amplitude  hatlos 

The  inverse  relaxation  length  for  the  thermal -neutron 
flux  in  the  assembly,  %   ,  wa3  determined  by  amplitude  ratios 
for  various  frequencies  of  oscillation.  The  value  actually 
obtained  was  the  inverse  attenuation  length,  but  as  previously 
noted,  under  certain  conditions  these  values  are  approximately 
equal.  The  experimental  conditions  were  such  that  the  above 
statement  is  true,  and  will  be  demonstrated  in  the  Discussion 
of  Results. 

Uncorrected -amplitude  ratios  were  used  in  place  of  foil 

-  X  z 
activities  in  the  equation  A  ■  Ce    .  With  the  method  of 

least  squares  an  initial  uncorrected  value  for  y  was  found 

(2,  p.  166).  Both  end -correction  and  harmonic -correct ion 

factors   were  determined  and  applied  to  the  amplitude  ratios 


* End-correction  factor:  Ce  -  l-e~2  * a{Q~z) 
Harmonic -correction  factor:  Ch»l-K^lie  " 
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"  *l3Z  +       *         e~  **iz 


%  IS  3^3  1 


2 mm  *aa£ 


*>33  J 

where     rmM  «  (*/&)*  (ma  +  n*  -  2)   +     Ylxe    . 
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with  this  Initial  value  of  y  .  The  origin  of  the  Z-axia  was 
taken  at  the  base  of  the  assembly}  therefore  hole  one  was  at 
a  position  of  z  •  12  inches.  A  second  value  for  y  was 
determined  with  the  corrected  amplitude  ratios.  This  process 
was  repeated  until  the  value  of  if   used  in  the  correction 
factors  agreed  with  the  If   found  by  using  the  corrected 
amplitude  ratios.  Computations  showed  that  harmonic  cor- 
rections were  significant  only  for  the  three  lower  holes  and 
end  corrections  for  the  two  upper. 

It  was  believed  that  harmonic  corrections  were  unnecessary $ 
therefore  a  value  for  if   was  determined  in  a  similar  manner 
using  the  end  corrections  only. 

D.  Determination  of  If   by  Median-flux  Levels 

Inverse  relaxation  length  was  also  determined  by  using 
the  median-flux  level.  Uncorrected  median -flux  levels  were 
used  in  place  of  foil  activities  in  the  equation  A  m  Ce"   2. 
A  procedure  was  followed  as  explained  in  the  preceding 
section  to  obtain  a  value  for  If   .  Values  were  determined 
by  us Ins  only  the  end -correction  factors. 

E.  Determination  of  Amplitude  Ratio 

Basically ,  the  need  of  the  amplitude  ratio  was  to  deter- 
mine if  the  amplitude  of  the  neutron  sine  wave  varied  with 
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the  frequency  of  the  source  oscillation.  To  eliminate  any 
"built-in"  variation  due  to  the  equipment,  the  calibration 
runs  next  to  the  oscillator  unit  were  used  as  the  basis  for 
determining  the  amplitude  ratios.  'Otoe   amplitude  ratio  was 
defined  as  the  amplitude  of  the  wave  obtained  for  a  particu- 
lar position  and  speed  divided  by  the  calibration-run 
aaplitude  at  the  same  speed.  Both  values  of  amplitude  were 
obtained  from  the  plots  of  the  data,  and  to  standardise  the 
results,  all  amplitudes  were  placed  on  a  one-minute  basis. 

After  the  amplitude  ratios  were  obtained  from  the  data, 
an  end-correction  factor  was  applied.  From  the  iterative 
process  explained  previously,  a  value  for  the  end -correct ion 
factor  was  obtained.  No  end-correction  factors  were  applied 
to  the  calibration  runs.  The  end-correction  factors  were 
divided  into  the  amplitudes  obtained  in  the  hole  positions, 
and  it  was  found  that  factors  were  significant  only  in  the 
top  two  holes. 

There  is  a  definite  difference  between  the  theoretical 
amplitude  ratio  obtained  from  Equation  31  and  the  amplitude 
ratio  obtained  from  the  experiment.  Theoretical  amplitude 
ratio  was  defined  as  the  ratio  of  amplitude  of  the  observed 
run  divided  by  that  of  an  infinitely  slow  run,  whereas  the 
experimentally-determined  amplitude  ratio  is  the  amplitude 
of  the  observed  run  divided  by  a  calibration  run.  Therefore, 
numerical  values  of  the  two  ratios  will  not  agree,  but  they 
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will  agree  in  showing  any  change  In  amplitude  ratio  with  fre- 
quency. 

F.  determination  of  Median-flux  Level 

The  median-flux  level  la  the  average  flux  at  a  particu- 
lar hole  position  and  oscillator  frequency.  From  the  plots 
of  original  data,  the  amplitude  and  the  lowest  value  of 
count  rate  found  in  the  graph  were  obtained.  The  raedian-flux 
level  was  then  taken  as  half  the  amplitude  added  to  the 
lowest  count  rate.  End -correction  factors  were  applied  to 
the  amplitudes  for  the  top  two  holes. 
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VII.  RESULTS 

The  original  data  for  all  runs  completed  including  the 
calibration  la  given  in  Tables  5  through  13  In  the  Appendix* 
Figure  5  shows  a  representative  plot  of  original  data  at  a 
frequency  of  2.31  radians/second  for  all  hole  positions  and 
calibration. 

Phase  angle  was  determined  for  each  run  and  the  results 
are  as  given  in  Table  1.  All  the  calibration  runs  proved 
that  there  is  no  "built  in"  phase  angle  in  the  equipment.  No 
phase  angle  is  observed  in  any  hole  position  above  0.63 
radians/second.  Approximately  a  20-degree  phase  angle  Is 
found  to  exist  in  all  five  hole  positions  at  0.63  radians/ 
second.  For  the  only  hole  position  investigated,  number  two, 
a  30-degree  phase  angle  Is  found  at  0.42  radians/second. 
Verification  of  the  phase  angles  in  hole  two  up  through  4.19 
radians/second  was  obtained  by  use  of  the  10-degree  cam.  Ho 
change  of  phase  angle  is  observed  in  proceeding  to  higher 
positions  while  taking  measurements  at  the  same  oscillator 
frequency. 

Amplitude  ratios  were  obtained  and  the  results  are  as 
given  In  Table  2.  Figure  6  is  a  plot  of  measured  amplitude 
ratio  versus  oscillator  frequency  for  hole  positions  one, 
three  and  five.  Amplitude  ratios  do  not  show  a  definite 
trend  with  frequency,  but  do  remain  approximately  constant 
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Table  1.  Phase  angle  (degrees) 


Frequency 
radians/second 

Hole  position 

Calibra- 
tion 

One 

Two 

Twoa 

Three 

Four 

Five 

0.42 

«••> 

30 

30 

0 

0.63 

20 

20 

20 

20 

20 

20 

0 

1. 

0 

0 

0 

0 

0 

0 

0 

1. 

0 

c 

0 

0 

0 

0 

0 

2.51 

0 

0 

0 

0 

0 

0 

0 

4.19 

0 

0 

0 

0 

0 

0 

0 

6.28 

0 

0 

•»»• 

0 

mm 

mm  mm 

0 

7.33 

0 

0 

MM* 

4MB 

mm 

vmmm 

0 

j.%2 

0 

0 

mm 

0 

12, 

0 

0 

mm 

0 

18. 

0 

0 

aRuns  accomplished  with  10-ogree  cam. 


Table  2.  Aaplitude  ratio 


Hole  position 

Frequency 

radians/second 

One 

Two 

Three 

Four 

Five 

0.63 

0.0710 
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Figure  6.  Amplitude  ratio  vs.  frequency 
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at  each  hole  position.  A  high  coefficient  of  variation  is 
evident  for  values  above  hole  one. 

Amplitude  ratios  were  also  plotted  against  vertical 
distance  in  the  assembly.  The  five  frequencies  investigated 
in  all  five  hole  positions  are  represented  on  this  plot, 
Figure  7»  The  line  slopes  drawn  correspond  to  the  y  values 
obtained  by  the  method  of  least  squares  for  the  set  of 
points  used. 

A  median-flux  level  for  each  hole  position  and  frequency 
n«s  determined.  Table  3  gives  the  results  of  this  computa- 
tion, which  also  includes  the  calibration  runs.  To  demon- 
strate the  change  in  median-flux  level  with  oscillator  fre- 
quency a  plot  of  these  two  variables  was  constructed  for  the 
first  three  hole  positions  and  calibration  run.  Figure  8 
shows  that  the  slopes  on  logarithmic  paper  of  these  lines 
for  the  four  positions  are  equal. 

Table  4  gives  the  various  values  of  y  obtained  by  using 
the  amplitude  ratios.  Uncorrected  values  gave  an  average 
of  0.0710  inches'"1  j  end-and -harmonic -corrected  values  gave 
an  average  of  0.0631  inches" %  and  end-corrected  values  gave 
an  average  of  0.0702  inches'-1.  End-corrected  2T  values  were 
plotted  versus  frequency  in  Figure  9*  Values  for  r   are 
approximately  constant  with  frequency.  Large  coefficients 
of  variation  are  noted  for  all  the  r   values  plotted  in 
Figure  9. 
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Table  3.  Median-flux  level  (total  counts) 


Frequency 
radians/ 

Hole  position 

Calibra- 

One 

fm 

Three 

Four 

Five 

tion 

second 

(2  min.) 

(2  mln. 

)  (3  mln,) 

(4  min.) 

(5  min.) 

(1  min.) 

0.63 
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3,535 

2,105 
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13,565 

1.03 

11,275 
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2,1 
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1. 

11,330 
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3,^73 

x,s 
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13,245 

2.51 
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1,098 
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mm 
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mm 
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Table  4.  Inverse-relaxation  length  (inches) 
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Figure  7.  Amplitude  ratio  vs#  vertical  distance 


kk 


00=  4.19 

RADIANS/SECOND 


VERTICAL    DISTANCE-INCHES 


^5 


20,000r 


CALIBRATION 


LU 

Z> 


^iopoo 

V- 

z 

O  7,000 

o 


LU 
> 

LU 


-HOLE    ONE 


4,000 


X 

Z> 
_J 

I 


Q 
LU 


HOLE     TWO 


2,000 


HOLE    THREE 


1,000 


0.3    0.4 


0.7         I  2  4  7         10 

FREQUENCY  -RADIANS/SECOND 


Figure   8.     Median-flux  level  vs.  frequency 


46 


io10 

LU 

X 

o 

E 

-p 

'    0.8 

— 

h- 

< 

> 

o 

t 

i 

i 

i 

z 

i 

i 

UJ 

_J 

i 

> 

0.6 

- 

z: 

g 

< 

-*- 

<0.4 

- 

_l 

LU 

or 

LU 

CO  0.2 

LU 

> 

z 

nn 

i       i      i     i    i    i                           i                i           ii 

0.4 


0.6         0.8        I  2  3 

FREQUENCY-  RADIANS/SECOND 


Figure  9.   Inverse  relaxation  length  vs.  frequency 


47 

A  value  for  tf   of  0,0661  inches*1  was  obtained  by  using 
the  method  of  median-flux  levels* 


4a 


VIII.  DISCUSSION  OF  RESULTS 

It  was  the  original  belief  that  a  phase  angle  would  be 
observed  at  any  oscillator  frequency  above  0*6  radians/second 
(Hi  p.  33)*  After  a  series  of  runs  had  been  completed,  this 
was  found  not  to  be  true.  Previously,  the  values  of  k  and  J? 
had  not  been  known  for  the  assembly  used,  and  had  to  be 
assumed.  This  led  to  the  erroneous  prediction. 

The  theoretical  analysis  was  performed  using  one  average 
group  of  delayed  neutrons.  Figures  10  and  11  show  the 
theoretical  phase  angle  and  amplitude  ratio  respectively  as 
determined  from  Equation  31*  A  value  for  k  of  0,53  was  used. 
With  f  equal  to  0*8?  (9,  p#  56),  and  Ba  equal  to  70  x 
10**  cms  (9#  P«  57b)  a  value  of  $.    was  determined  to  be 
C. 00090  seconds  from  the  equations 

L*  -  ytW )  and  '  «  x  *fhV 

and  an  assumed  value  of  i0  •  0.001  seconds.  The  total 
fraction  (3  of  the  delayed  neutrons  was  set  equal  to  a* 0064 
(7)  with  a  corresponding  value  of  X    equal  to  0.08  seconds*1. 
It  can  readily  be  n«mn  from  Figure  10  that  a  phase  angle 
would  not  be  observed  in  the  assembly  until  an  oscillator 
frequency  of  40  radians/second  is  reached.  A  phase  angle 


*Beck,  0.  It* *  Ames,  Iowa.  Data  from  partially  completed 
M.  s.  Thesis.  Private  cotesunication.  1959. 
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does  not  appear  for  this  investigation  until  very  high  fre- 
quencies because  the  value  of  k  is  low.  If  k  were  0*9  a 
phase  angle  would  be  observed  at  a  frequency  of  20  radians/ 
second*  The  overall  shape  of  the  theoretical  curve  regains 
the  same  with  varying  values  of  k. 

Huns  were  completed  to  the  lower  and  upper  frequency 
limits  of  the  oscillator  and  counting  equipment.  The  upper 
limit  (19  radians/second)  is  short  of  the  required  frequency 
for  an  observable  phase  angle.  In  complete  agreement  with 
theory,  no  phase  angle  was  observed  in  any  of  the  hole  posi- 
tions above  Q.63  radians/second.  Use  of  the  10 -degree  cam 
verified  these  results. 

Readings  were  taken  in  all  hole  positions  for  a  frequency 
of  0*63  radians/second.  A  consistent  phase  angle  of  approxi- 
mately 20  degrees  is  found.  At  0.42  radians/second  the  phase 
angle  in  hole  two  is  approximately  30  degrees.  The  theory 
advanced  does  not  predict  a  phase  angle  at  these  low  values 
of  frequency.  Further  investigations  were  not  conducted  into 
this  low  frequency  range  since  the  equipment  did  not  provide 
adequate  speed  control  below  4  epm. 

Theory  makes  no  distinction  of  position  in  relation  to 
phase  angle |  therefore  the  phase  angle  should  be  identical 
in  all  areas  of  the  assembly,  Ko  change  in  phase  angle  with 
hole  position  was  observed  during  the  investigation.  Surveys 
were  taken  only  up  the  center  of  the  assembly  so  the 
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proof  of  the  constancy  of  phase  angle  la  limited  to  this  con- 
dition. With  the  method  of  source  forcing  used  It  is  believed 
that  some  large  variations  would  be  found  at  the  edge  of  the 
assembly  near  the  base.  Simply  noting  that  the  pattern 
position  facing  these  extreme  positions  will  vary  as  much 
as  45  degrees  from  the  recorded  position  indicates  the 
physical  reason  for  the  nonconformance. 

As  shown  in  Figure  5*  the  curve  of  the  neutron  flux 
measured  at  the  various  hole  positions  is  shown  to  be 
sinusoidal.  Scatter  of  points  increases  as  higher  hole 
positions  are  investigated.  *Jhe  range  of  frequency  for 
usable  data  narrows  as  measurements  are  taken  higher  In  the 
assembly.  Statistical  variations,  and  inaccuracies  in  the 
detector  and  scaler  distort  the  counts  received.  Due  to  the 
already  small  amplitude  of  the  sine  wave,  any  distortions 
cause  a  scattering  of  points  that  make  the  readings  unusable. 
This  situation  Is  aggravated  by  higher  frequencies.  As  the 
oscillator  frequency  increases,  the  neutron  waves  are  propa- 
gated quite  close  together.  Any  disturbance  or 
irregularity  in  the  entire  assembly  or  counting  system  will 
cause  an  incorrect  measurement  to  be  taken.  For  example,  a 
measurement  is  to  be  taken  at  a  specified  oscillator-unit 
degree.  A  small  disturbance  or  Irregularity  will  cause  the 
neutron  wave  to  be  delayed  in  reaching  the  detector.  At  a 
high  frequency  this  delay  is  significant  since  the  portion 
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of  the  sine  wave  actually  surveyed  may  be  30  to  40  degrees 
removed  from  the  portion  desired* 

As  shown  In  Figure  5»   the  standard  deviation  at  hole  one 
is  reasonable  and  acceptable.  However  as  readings  are  taken 
at  higher  positions  in  the  assembly  the  coefficient  of  varia- 
tion caused  by  standard  deviation  increases  even  though 
longer  counting  times  are  allowed.  This  is  caused  by  lower 
count  rates  due  to  the  exponential  attenuation  of  the  thermal - 
neutron  flux.  An  average  coefficient  of  variation  for  all 
the  runs  completed  ranges  from  one  to  three  percent  In  going 
from  hole  one  to  five.  It  is  seen  then  that  to  obtain 
eoeff icients  of  variation  of  less  than  two  percent  at  the 
higher  hole  positions,  very  long  counting  times  are  needed, 
in  the  order  of  10  minutes  at  hole  five. 

Amplitude  ratios  were  determined  and  are  given  in  Table 
2.  These  amplitude  ratios  for  holes  one,  three  and  five  were 
plotted  against  frequency  in  Figure  6.  A  least  squares  analysis 
was  applied  to  the  data  for  holes  three  and  five,  and  it  was 
found  that  the  line  joining  hole-three  data  had  a  slope  of 
-0.1d4  *0.22$,  and  the  line  joining  hole-five  data  had  a  slope 
of  -Q.231  20.34*f.  This  Indicates  that  a  horizontal  line 
representing  a  constant  value  of  amplitude  ratio  with  frequency 
lies  easily  within  the  standard  deviation  for  the  slopes  found 
for  both  holes.  It  is  believed  that  the  amplitude -ratio  values 
are  approximately  constant  as  predicted  by  Figure  11  and  that 
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the  slopes  shown  on  Figure  6  are  caused  by  the  high  standard 
deviation  rather  than  variance  of  amplitude  ratio  with  fre- 
quency. 

Figure  11  has  the  overall  shape  that  the  amplitude  ratio 
versus  frequency  would  have  for  any  value  of  te.  For  a  k 
greater  than  0,53,  an  amplitude  ratio  less  than  one  would  be 
noted  at  frequencies  less  than  indicated  on  Figure  11. 

Scatter  of  the  points  in  Figure  6  is  due  to  a  very  high 
coefficient  of  variation.  As  was  mentioned  previously*  some 
scatter  of  points  is  present  at  all  positions  especially  in 
the  higher  hole  positions  at  higher  frequencies.  Error  occurs 
when  fairing  in  the  curves  to  fit  the  data.  A  very  signifi- 
cant error  occurs  when  the  standard  deviation  of  the  points 
is  large  compared  to  the  amplitude*  For  example,  at  hole 
three  a  standard  deviation  of  approximately  60  counts  is 
foundj  however  an  amplitude  of  only  300  counts  exists. 
Coefficient  of  variation  caused  by  standard  deviation  of  the 
amplitude  is  then  approxliaately  20  percent.  The  condition  is 
exaggerated  at  higher  hole  positions  such  that  a  36  to  42 
percent  coefficient  of  variation  in  hole  five  is  noted. 
Significant  errors  appear  for  the  amplitude  ratios  above 
hole  one.  Considerably  longer  counting  times  are  needed  at 
all  higher  hole  positions  to  give  satisfactory  results. 
A  counting  time  of  approximately  10  minutes  would  give  a  13 
percent  coefficient  of  variation  for  hole  four. 
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Bven  though  the  coefficients  of  variation  for  the  ampli- 
tude ratios  are  quite  high,  the  results  are  believed  to  be 
sufficiently  accurate  and  repetitive  to  support  the  conclu- 
sions drawn. 

Values  of  the  medlan-f lux  level  are  given  in  Table  3  and 
the  plot  of  level  versus  frequency  is  shown  as  Figure  8. 
Deviation  was  so  slight  in  holes  four  and  five  that  they  were 
not  included  in  the  figure. 

The  transfer-function  development  does  not  predict  a 
decrease  in  the  median  flux,  yet  all  four  curves  show  that 
median-flux  level  does  decline  as  frequency  increases.  Slopes 
on  the  logarithmic  plot  (Figure  6)  of  all  four  lines  are 
found  to  be  equal.  If  each  line  is  represented  by  the  equa- 
tion 0  »  K  CO  n,  where  K  is  a  constant  and  n  is  the  slope  of 
the  line,  then  the  rate  of  neutron-flux  decline  with  fre- 
quency is 

d#     Jin  co  n     n# 
&u  to  cj 

The  relative  rate  of  decline  is 

a/     n 

&CJi     m  CO 

9 
With  identical  values  of  n,  and  the  same  values  of  to ,  each 
hole  position  has  the  same  relative  rate  of  decline.  The 
decline  is  therefore  not  a  function  of  the  assembly. 
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The  decline  is  attributed  to  either  the  method  of  source 
forcing  or  to  the  counting  equipment.  No  logical  reason  can 
be  advanced  for  the  unit  itself  to  cause  the  decline ,  but 
the  counting  equipment  could  easily  have  caused  this  un~ 
predicted  result,  The  number  of  activations  of  the  preset- 
count  relay  is  in  proportion  to  the  frequency.  A  small  relay 
lag  would  then  become  important  at  higher  frequencies  and 
significantly  reduce  the  amount  of  time  the  scaler  was  count- 
ing* Furthermore ,  investigation  shows  that  the  reduction  of 
thermal  neutrons  and  fast  neutrons  counted  is  In  proportion 
to  their  total  counts.  The  deviation  from  theory  is  there* 
fore  said  to  be  a  mechanical  irregularity  of  the  counting 
equipment  rather  than  a  function  of  the  assembly. 

Amplitude  ratios  were  not  affected  by  the  decline  of 
counts  with  frequency  since  the  calibration  runs  were  used 
as  the  denominator  of  the  ratio. 

When  plotted  on  semi-log  paper,  amplitude  ratio  versus 
vertical  distance  gives  a  reasonably  straight  line.  Figure  7 
shows  the  results  for  frequencies  from  0.63  to  4.19  radians/ 
second.  As  noted  previously,  the  large  coefficients  of 
variation  of  data  for  the  higher  hole  positions  result  In 
some  scatter  of  the  amplitude  ratios.  It  is  seen,  however, 
that  the  points  nearly  all  lie  close  to  the  straight  lines 
indicating  that  the  points  can  be  considered  reliable. 
Standard  deviation  of  the  amplitude  ratios  is  indicated  on 
the  figure. 
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Values  for  inverse  relaxation  length,  tf"  ,  were  computed 
for  the  frequencies.  These  values  obtained  by  an  iterative 
process  are  listed  in  Table  4,  and  it  is  found  that  values 
range  from  0.0i>82  to  0.0729  inches* x.     A  value  for  X   was 
found  without  using  either  the  end  or  harmonic  corrections; 
the  average  value  was  0.0710  inches*1.  When  both  corrections 
wei*e  applied,  an  average  value  of  0.0631  inches"1  was  found. 
It  is  believed  however  that  the  higher  harmonics  have  been 
greatly  iesseneu  with  this  particular  arrangement  of  sources 
in  the  oscillator.  The  sources  are  placed  in  a  line  and  are 
surrounded  hy   three  inches  of  paraffin  which  thermal izes 
many  of  the  neutrons  produced.  The  paraffin,  lined  sources, 
and  two  water-filled  tanks  in  the  outermost  pedestal  spaces 
help  to  scatter  the  neutrons  more  uniformly  beneath  the  base, 
so  chat  a  plane  source  of  thermal  neutrons  is  approximated. 
Dopchie,  Leonard,  Neve  de  Mevergnies  and  Tavernler  (3) 
showed  in  their  work  with  an  exponential  assembly  that  dis- 
tributing four  sources  in  a  square  form  reduced  the  harmonic - 
correction  term  to  |  percent  or  one-fourth  the  correction 
needed  when  the  sources  were  placed  together. 

With  application  of  just  the  end  corrections  an  average 
value  for  2T  of  0.0702  inches"1  was  determined.  Hayes  (o, 
p.  27)  found  by  the  method  of  foil  irradiation  that  the 
assembly  in  this  configuration  had  a  2f  of  0.0705  inches*1. 
Harmonic -correction  factors  would  appear  then  not  to  be 


needed,  or  if  used  the  full  correction  should  not  be  applied. 
Line  slopes  on  Figure  7  are  the  X  values  determined  by  using 
only  the  end  corrections* 

Values  of  inverse  relaxation  length  obtained  by  using 
end  corrections  only  were  plotted  with  frequency  on  Figure  9* 
A  standard  deviation  for  V   was  obtained  by  using  the  method 
of  propagation  of  precision  indexes  as  presented  by  worthing 
and  Geffner  (1?,  p.  209)  and  are  indicated  on  the  figure. 
Large  coefficients  of  variation  are  again  noted  for  2f   • 
This  is  due  to  the  high  coefficients  of  variation  of  the 
amplitude  ratios,  especially  at  the  higher  hole  positions* 
By  using  longer  counting  times  the  relative  standard  devia- 
tion of  the  amplitude  ratios  would  be  greatly  reduced  and 
statistically  better  values  of  X   would  follow*  For  example, 
if  an  average  relative  standard  deviation  of  5  percent  were 
found  for  the  amplitude  ratios*  then  the  coefficient  of 
variation  for  ?   in  this  experiment  would  decrease  by  a 
factor  of  approximately  six* 

It  is  seen  from  Figure  9  that  there  is  no  trend  in  y 
with  increasing  frequency*  A  constant  value  of  V   should  be 
found,  but  due  to  the  high  standard  deviations  values  differ- 
ing slightly  from  the  average  of  0.0702  inches* l   are  noted* 

To  substantiate  the  values  of  X   determined  by  using 
amplitude  ratios  an  average  value  of  X   was  found  by  using 
the  median-flux  levels.  As  shown  on  Figure  12  the  values  of 
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the  median  flux  are  grouped  close  together*  An  average  value 
was  obtained  from  each  group  ana  a  value  of  X   computed  from 
them  by  using  end  corrections  only.  Standard  deviation  of 
the  average  points  is  approximately  the  else  of  the  symbols. 
therefore  the  value  of  ¥   obtained  is  reliable,  rrora  Figure 
12  a  value  of  X   was  found  to  be  COool  inches"1  in  compari- 
son to  the  average  value  of  0.0702  inches'"1  computed  by  using 
amplitude  ratios.  Agreement  within  3  percent  between  the 
two  methods  indicates  the  amplitude -ratio  method  to  be 
acceptable • 

The  validity  of  the  statements  concerning  the  attenua- 
tion length  and  l/y  was  investigated.  Values  used  in  the 
computation  were  T  «  350  cm*  (16,  p.  331)#  kp  •  k(l  -  (3  )  « 
0.526,  L3  m  325  cm*,  i0  •  0.001  seconds  and  CO  m  20 
radians/second.  A  value  was  determined  from  Equation  32  for 
attenuation  length  and  found  to  be  14.8  inches.  The  inverse 
of  this  number  is  0.0676  inches"1  which  agrees  favorably 
with  the  values  of  0.0702  inches*4  and  0.06^1  inches"1  ob- 
tained by  the  methods  oS  amplitude  ratios  and  median-flux 
levels  respectively. 

Frequencies  utilised  in  the  experiment  were  small  enough 
to  cause  the  attenuation  length  to  be  approximately  equal 
to  the  relaxation  length.  From  Equation  St  it  is  noted  that 
a  frequency  of  20  radians/second  produces  ne&iiglbe  change 
in  the  value  of  attenuation  length  as  compared  to  very  low 
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frequencies.  Hence  from  the  previous  computation*  and  dis- 
cussion it  appears  that  the  attenuation  length  and  conse- 
quently the  complex  buckling  for  the  frequencies  investigated 
are  approximately  the  same  as  relaxation  length  and  material 
buckling*  Therefore  it  would  be  proper   in  this  case  to 
utilise  the  values  of  attenuation  length  for  determining  y  . 
Even  with  the  combination  of  varying  and  steady  neutron  flux, 
computation  of  ^  by  amplitude  ratio  and  median -flux  level 
is  feasible  since  the  frequency  n&vev   is  large  enough  to 
produce  significant  changes  in  X   » 

For  assemblies  with  much  larger  values  of  kp,  frequencies 
used  in  this  experiment  would  produce  changes  in  the  attenua- 
tion length  and  give  varying  values  for  the  complex  buckling, 
therefore  the  inverse  attenuation  length  would  not  be  equal 
to  the  inverse  relaxation  length. 
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IX.   CONCLUSIONS 

A  theoretical  and  experimental  investigation  into  neutron 
diffusion  in  the  subcritical  assembly  with  sinusoidal  source 
oscillation  was  conducted  at  Iowa  State  College*  An  oscilla- 
tor unit  produced  the  sinusoidal  source  oscillation  of  the 
thermal  neutrons,  while  a  selective -count  device  allowed 
discrete  sections  of  the  sine  wave  to  be  observed  and  counted. 
Neutron-diffusion  information  was  extracted  from  plots  of 
these  counts.  Conclusions  drawn  from  the  information  are  that; 

1.  No  phase  angle  is  observed  for  oscillator  frequencies 
between  O.63  and  !$«&$  radians/second.  This  is  in  agreement 
with  theory. 

2.  An  unpredicted  phase  angle  of  20  degrees  is  observed 
at  a  frequency  of  0.42  radians/second  and  an  angle  of  30 
degrees  at  0.63  radians/second. 

3.  There  Is  no  phase-angle  change  at  higher  hole  posi- 
tions on  the  centerline  of  the  assembly  at  the  same  oscillator 
frequency,  which  ranged  from  O.63  radians/second  to  6.28 
radians/second.  This  is  in  agreement  with  theory. 

4.  Amplitude  ratios  for  a  given  hole  are  approximately 
constant  for  varying  frequencies  from  O.63  radians/second 

to  6.28  radians/second.  This  is  in  agreement  with  theory. 

5.  The  median-flux  level  decreases  as  frequency  in- 
creases. Mechanical  operation  of  the  counting  equipment 
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causes  this  change  rather  than  any  assembly  peculiarities  or 
characteristics.  A  change  in  level  is  not  predicted  by 
theory* 

6*  She  oscillator  unit  thermallzes  enough  neutrons  to 
give  satisfactory  readings  up  through  hole  five  up  to  4*19 
radians/second.  Any  frequency  readings  above  12*35  radians/ 
second  are  confined  to  hole  one. 

7*  An  inverse  relaxation  length  of  0*0702  inches"1 
nee  determined  by  using  amplitude  ratios. 

§•  An  inverse  relaxation  length  of  0.0651  inches*1  was 
determined  by  using  median-flux  levels. 

9.  Tlie  attenuation  length  of  the  neutron  wave  is  equal 
to  a  constant  value  of  l/'y  for  the  assembly  and  frequencies 
investigated. 

10*  Longer  count  times  are  needed  to  reduce  coefficients 
of  variation  in  the  amplitude  ratios*  In  all  other  respects 
the  experimental  equipment  and  procedure  will  yield  satis* 
factory  results  within  certain  frequency  limits  for  phase 
angle,  amplitude  ratio  and  inverse  relaxation  length. 
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X.  SU30ESfI0NS  FOR  FOHTIffiR  STUB* 

From  a  review  of  the  literature  it  is  obvious  that  neutron 
diffusion  in  a  subcrltical  assembly  with  sinusoidal  source 
oscillation  has  not  been  fully  investigated. 

several  ways.  Another  speed  reducer  could  be  introduced  to 
allow  speeds  below  4  rph  to  be  utilised.  2he  cam  could  be 
replaced  by  a  phototube  device  that  would  trigger  the  scaler 
at  the  higher  frequencies.  An  automatic  revolution  counter 
would  be  necessary  to  measure  the  frequency  at  the  higher 
speeds.  A  gating  circuit  would  be  more  desirable  to  inter- 
rupt pulses  going  to  the  scaler  than  using  the  preset-time 
relay  which  is  burdened  at  high  frequencies.  Still  another 
possibility  might  be  to  use  an  ionization  chamber  connected 
to  a  sensitive  electrometer.  The  wave  produced  by  a  re- 
corder connected  to  the  electrometer  could  then  be  analysed 
for  data  desired. 

In  this  experiment,  only  the  center  of  the  assembly  was 
Investigated.  Measurements  could  be  taken  throughout  the 
assembly  to  see  if  the  thermal-neutron  flux  behaved  in  the 
same  manner  as  the  flux  produced  by  a  steady  source.  Further 
theoretical  and  experimental  work  with  neutron  waves  could  be 
accomplished  with  an  investigation  into  the  change  in  prompt 
buckling  by  source  forcing.  Speeds  considerably  below  and 
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above  those  Investigated  here  could  also  be  utilised  to  ex- 
pand the  experimental  verification  of  the  theoretical  transfer 
function,  the  above  suggestions  could  be  expanded  greatly 
by  changing  the  lattice  In  the  assembly  and  repeating  the 
particular  Investigation.  Voids  and  absorbers  could  also 
be  placed  In  the  assembly  to  observe  the  effect  of  them  upon 
the  neutron  wave. 


XI.   LITERATURE  CITED 


1.  Axtmann,  K.  C,  Dessaucr,  3.  and  Parkinson,  T.  P. 

Reactivity  measurements  in  a  suberitieal  pile. 
U.  S.  Atomic  Energy  Commission  Report  DF-48 
(Du  Pont  de  Nemours  (E.  I.)  and  Co.,  Wilmington, 
De 1 . ]  November  1955* 

2.  Dana,  ?.  C.  and  Killyard,  L.  R.  Engineering  problems 

manual.  New  York,  N.  Y.,  McGraw-Hill  Book  Co., 
Inc.  19^7. 

3.  Dopchie,  K.,  Leonard,  P.,  Have  de  Mevergnies,  M.  and 

Tavernier,  G.  Conducting  an  exponential  experiment 
with  a  natural  -  U  graphite  lattice,  Nucleonics 
Ihi   37-60.  March  1936. 

4.  Franz,  J.  P.  Pile  transfer  functions.  U.  S.  Atomic 

Energy  Commission  Report  AECD-3260.   [Technical 
Information  Service,  AEC].  July  13,  19^9. 

5.  Glasstone,  S.  Principles  of  nuclear  reactor  engineer- 

ing. Princeton,  N.  J.,  S«  Van  No© trend  Co.,  Inc. 
1955. 

6.  and  Edlund,  n.  C.  The  elements  of  nuclear  re- 

actor  theory.  Princeton,  N.  J.,  D.  Van  Nostrand 
Co.,  Inc.  1952. 

7.  Hadley,  J.  W.  Speech  delivered  to  the  Iowa  State  College 

American  Nuclear  Society  Student  Chapter.  Univ.  of 
Calif.,  Livermore,  Calif.,  Author.  January  15, 
1959. 

8.  Mayes,  J,  T.  Flux  distribution  in  a  unit  cell  of  a 

uranium  graphite  suberitieal  assembly.  Unpublished 
M.S.  Thesis.  Ames,  Iowa,  Iowa  State  College 
Library.  1956. 

9.  fioganson,  J.  H.  Operating  characteristics  of  a  uranium 

graphite  suberitieal  assembly  with  coolant  simula- 
tion. M.S.  Thesis.  Ames,  Iowa,  Iowa  State  College 
Library.  1957* 

10.  Martin,  D.  H.  Correction  factors  for  measurements  with 
cadmium  covered  foils.  U.  S.  Atomic  Energy  Com- 
mission Heport  NAA-SR-1076  [Hcrth  American  Aviation, 
Inc.,  Downey,  Calif.].  October  15,  195^. 


67 

11*  Murray,  K.  L.  Nuclear  reactor  physics* 

Cliff a.  N.  J.,  Prentice-Hull,  Inc.  1957 

12.  Owens,  J.  I.,  Jr.,  Crever,  F.  E.  and  Pigott,  J.  H. 

A  proposed  automatic  control  system  for  the 
Schenectady  reactor.  U.  S.  Atomic  Snergy  Coasaission 
Report  ASOD-4209  [Technical  Information  Service , 
A£C]«  $areh  29*  19*9. 

13.  Ricci,  «,  J«  Transfer  function  of  a  uraniua-graphite 

subcritlcal  assembly.  Unpublished  M.S.  Thesis. 

Ames,  Iowa,  Iowa  State  College  Library.  1. .. 

14.  Soodak,  H.  Pile  kinetics.  In  woodman,  C,  ed.  The 

science  and  engineering  of  nuclear  power.  Vol.  2, 
pp.  09-102.  Cambridge,  Mass.,  Addiaon-Wealey 

Press  *  19^9 • 

15.  Weinberg,  A.  N*  and  Schweinler,  K.  C.  Theory  of 

oscillating  absorber  in  a  chain  reactor.  Physical 
Review  7**»  3^1*863.  1 

**«    and  Winner,  E.  P.  The  physical  theory  of 

neutron  chain  reactors.  Chicago,  HI.,  The  Univer- 
sity of  Chicago  Press.  191. 

17.  Worthing,  A.  8.  and  Oeffner,  J.  Treataaent  of  experi- 
mental data.  New  York,  H.  Y.,  Colm  Wiley  and  Son.*, 
Inc.  1 


The  author  wishes  to  thank  3r.  Robert  £.  Unrig  for  his 
original  suggestion  of  this  project}  and  especially  for  his 
tireless  assistance  and  sincere  interest  in  all  phases 
encountered  both  in  this  specific  project  and  in  the  course 
work  here  at  Iowa  State  College*  Many  thanks  are  due  Dr. 
Glenn  Murphy  for  his  guidance  and  encouragement  given  to 
the  author  during  this  course  of  study. 

This  thesis  culminates  three  years  of  study  under  the 
auspices  of  the  U.  S.  Navy  in  its  postgraduate  training 
program.  The  author  wishes  to  express  his  appreciation  to 
the  United  States  Naval  Postgraduate  School  for  making  this 
year  at  Iowa  State  College  possible. 


XIII.     APFENDXX 


70 


Table 


j» 


Original  data  at  0*63  radians/second  (total 

number  of  counts*) 


Oscillator 

] 

tola  position 

position 

One 

Tuo 

Three 

lour 

Five 

degrees 

(3  «ain.) 

(3  min.) 

(3  min.) 

(4  min.) 

(5  min.) 

0 

17,900 

3*7 

2,254 

1,234 

6o 

10,511 

3, 

2,220 

1,257 

16,310 

8,600 

3,o60 

2,170 

1,231 

9C 

l. 

L§*c 

2,105 

1,199 
1,134 

120 

17,276 

^,189 

3,6^0 

2,102 

130 

16,306 

7,847 

3,522 

2,031 

1,074 

loO 

16,c. 

7,708 

3,374 

I'M 

1,042 

210 

13,595 
16,2* 

7,528 

3,^ 

1,100 

240 

7,7^5 

3,484 

2,016 

1,077 

270 

17,0 

7,  pa 

3,491 

2,036 

1,086 

300 

16,971 

7,625 

3,496 

2,040 

1,103 

330 

17,747 

3o 

3,665 

2, 11C 

1,164 

Table  6.  Original  data  at  1.05  radians/second  (total 
number  of  counts) 


Oscillator 

position 

degrees 


(2  min.) 


Hole  position 

Two     Three     ?our     Five 
(2  min.)  (3  rain.)  (4  min.)  (5  sain.) 


g 

ISO 

ISO 
210 
»4C 
270 
300 
330 


11,950 

6,342 

3,963 

11,. 

3,937 

11, col 

5,721 

3,716 

11,471 

5,386 

3,715 

11,017 

5,351 

3,701 

10,7^9 

5,230 

3,498 

10,455 

5,223 

5,193 

3,5&2 

11,207 

5,233 

3,660 

11,777 

5,2b7 

3,554 

11,210 

3,769 

11,373 

5,7 

3, 

2,1 

2,165 
2,215 
2,193 

!:g| 

1,984 

2,034 
2,111 
2,206 

2,239 


1,297 
1,370 

1,309 

1,282 
1,201 

1,2-09 
1,200 
1,242 

1,268 
1,229 
1,249 
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Table  7.  Original  data  at  1.68  radians/second  (total 
number  of  counts ) 


Oscillator 
position 

Aegrees 

Hole  position 

One      Two     Wiree     Four     Five 
(2  rain.)  (2  rain.)  (3  min.)  (4  rain.)  (5  mln.) 

90 
lac 

210 
240 
270 

300 
330 


ia,i7S 

5,510 

3,563 

12,061 

5,513 

3,479 

11,0^2 

5,426 

3,44a 

11,471 

5,2o3 

3,540 

11,065 

5,282 

3,444 

10.59O 

£,114 

3,460 

10,62f> 

4,9^2 

3,310 

10,571 
10,94O 

5,000 

3,439 

5,214 

3,369 

11,291 

5,191 

3,525 
3,579 

11,665 

5,189 

12,257 

5,234 

3,499 

2,082 

2,090 
2,154 

2,067 

1,959 
1,932 

1,924 
1,394 

1,951 
1, 

1,956 
1,98a 


1,110 

1,151 

1 ,096 

1,099 
1,113 
l,r 
1,076 

1*035 

1,028 

1,039 

l,llu 

1,131 


fable  b.  Original  data  at  2*51  radians/second  (total 
nuraber  at   counts) 


Oscillator 

Hole  position 

position 

One 

Two 

Three 

Four 

Five 

degrees 

(2  rain.) 

(2  rain.) 

(3  rain.) 

(4  rain.) 

(5  min.) 

0 

12,027 

5,431 

3,656 

2,196 

1,235 

3° 
Sc 

11,807 

5,407 

3,641 

2,110 

1,221 

11,729 

5,42o 

3,588 

2,177 

1,212 

90 

11,224 

5,3^4 

3,5i>9 

2,( 

1,1 

120 

10,765 

5*381 

3,522 

2,034 

1,141 

150 

10,467 

5,197 

3#;- 

2,000 

1,151 

180 

10,315 

5,153 

3,394 

2,002 

1,110 

210 

10,37? 

10,814 

3,252 

3,453 

1,973 

1,124 

240 

5,2 

3,526 

2,092 

1,149 

270 

11,421 

5,333 

5,424 

3,801 

2,150 

1,1 

300 

11,515 

3,584 

2,136 

l,lb6 

330 

11,391 

5,436 

3,673 

2,194 

1,169 
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Table  9*  Original  data  at  4.19  radians/second  (total 
number  of  counts) 


Oscillator 

Hole  position 

position 

One 

Two 

Three 

Pour 

ive 

degrees 

( 2  min . ) 

(2  min.) 

( 2  min . ) 

(4  min.) 

(5  min.) 

0 

11,590 

5,316 

2,417 

2,073 

1,124 

00 

11,585 

5,406 

2,396 

2,053 

1,145 

11,376 

5,307 

2,319 

1,974 

1,142 

90 

10,395 

10,44d 

5,23t> 

2,337 

1,967 

1,126 

150 

5,326 

2,279 

1, 

1,063 

150 

10,209 

4,975 

2,221 

1,944 

1,072 

160 

10,137 

4,857 

2,183 

1, 

1,032 

210 

10,269 

4,696 

2,241 

1,949 

1,067 

240 

10,6; 

27 

2,252 

1,952 

1,127 

270 

11,423 
11,246 

4,992 

2,266 

2,001 

1,108 

300 

5,222 

2,309 

1,994 

1,113 

330 

11,360 

3,406 

2,3b5 

2,103 

l,12o 

Table  10, 

Original  data  at  6.28 

radians/second  (total 

number  of  counts) 

Oscillator 

Mole  position 

position 

One 

Two 

Three 

degrees 

(1  min.) 

(1  min, ) 

(3  min. ) 

0 

5,343 

3,211 

4,086 

30 

5,706 

3,241 

4,051 

So 

5,<- 

3,156 

4,012 

90 

mum 

3,12C 

3*S 

120 

5,641 

3,013 

3,o57 

150 
ISO 

5,337 

3,B23 

5,271 

2,9& 

320 

210 

5,334 

3*050 

3, 

240 

5,521 

3,129 

3,- 

270 

-«. 

3#l 

3*S97 

300 

5,720 

3, loo 

|»99C 

33o 

5,735 

J,24p 

4,030 
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Table  11.  Original  data  at  hole  petition  one  (total  number 
of  counts) 


Oscillator 

Radians/sec  ond 

position 

,  7-33  t 

,42 

12. _ 

18.55 

degrees 

(1  *&&•* 

(2  &In.) 

(1  mln.) 

(1  rain.) 

0 

16,3*7 

11,2a 

5,367 

3° 

16,47-. 

11 ,253 

5, 

5,oU 

oO 

15,9.. 

10,843 

5,414 

5,3^9 

90 

1%U 

1C3 

5, 

5,1^1 

120 

15,2- 

10 . 

5,0 

4,  6$4 

(35 

jO 
160 

13,1* 

Kb 

79 

13,2. 

10,140 

H.a 

4,700 

210 

15,3: 

*#s 

240 

15,932 

10,729 

5,1 

5,241 

270 

16 f 402 

10,970 

5,218 

5,326 

300 

16,3 

11,133 

5,508 

330 

16,655 

11,1 

5,  05 

5,; 

Table  12. 

Original  data  at  hole  position  two  (total  number 

of  counts) 

Oscillator 

position 

Hadlana/second 

7.33 

9.42a 

9.^2b 

12.55 

degrees 

{ 3  min , ) 

(2  mln.) 

( 2  rain . J 

(1  rain.) 

0 

7,51 

5,1 

5,1-2 

51 

30 

7,^ 

5,- 

5,197 

4,906 

60 

7,602 

$#< 

5,021 

4,97^ 

w 

7,3&2 

5,< 

120 

7,29c 

4,9 

5,110 

4,647 

150 

7,; 

^,975 

5,111 

4,634 

1 

§ , ."  ■ 

^,573 

21 

7,3 

*,- 

^,972 

*,751 

240 

7,: 

V975 

5,3 

4,7^0 

270 

7,67V 

5,1*7 

5,239 

X) 

7,331 

5,  J 

5#J 

^,938 

330 

7,539 

5,232 

5,270 

5,063 

aFirst  run. 
^Second  run. 
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Table  13.  Original  data  for  calibration  (counts/minute) 

Oscillator 

position  Radians/second 

degrees 0,42  O-63    1.03  1.66  2.51  4.11- 

0                   16,6^6  13,315     1  15,700  18,218  ic;,003 

30                  17,37  -    ,291     1  )  17,513 

16,174  17,270     17,366  17,475  17,377  17^ 

14,729  >0     16,140  1.  15.615  14, 

Ifl                     13,227  13,609     13,153  12,384  12,576  12,442 

IJO                    9,155  9,543           701  9.:  9,433  ?57 

lao                   3,049  7,9*7      h&°  T.T90  7,677  7,460 

10                    o,573  -,-                >51           -07  -,371 

40                  10,693  11,497     11,333  10,939  10,955  10,726 

270                    —  14,412     14,557  14,104  13,881  13,591 

300                    —  17,070    16,562  1   ,  16,343  16,223 

330                    —  18,3h2     18,198  17,777  17,763  17, 


6.;  7.  12. 

~— — — — — -  —     ■»■    ■■       '■  ■■'      — ■-  ■■-} ^-r-- — , , 1 1 ^-p-p— — 1 -  ■         , «.. ^-^~ ^  tfr.    .tfl    m|1   ..|f. 

0                  17,694  17,097  17,353  17,3*3  17,162 

30                  16,  17,701  16,959  17,12o  16, 

60                 15,703  16,3^1  15,450  15,012  15, 

90                   13,244  14,206  12,912  13.748  13,023 

120                     10,517  11,001  10,019  10,790  10,726 

150                   7,783  -,505  7**07  8,*  7,923 

7,171  7,.                >39  6,:  6,709 

210                       o,419  7,758             01  0,144 
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